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The QCD coupling a s is extracted from the energy evolution of the first two moments (mul¬ 
tiplicity and mean) of the parton-to-hadron fragmentation functions at low fractional hadron 
momentum z. Comparisons of the experimental e + e~ and deep-inelastic e ± p jet data to our 
NNLO*+NNLL predictions, allow us to obtain a s (m 2 ) = 0.1205±0.0010lg'gggg, in excellent 
agreement with the current world average determined using other methods at the same level 
of accuracy. 


1 Introduction 

For massless quarks and fixed number of colours, the only fundamental parameter of the theory 
of the strong interaction, quantum chromodynamics (QCD), is its coupling constant cr s . Starting 
from a value of A QCD ~ 0.2 GeV, where the perturbatively-defined coupling diverges, o s decreases 
with increasing energy Q following a 1/ln(Q 2 /A 2 CD ) dependence. The current uncertainty on 
o s evaluated at the Z mass, o s (m 2 ) = 0.1185±0.0006, is ±0.5% 4 , making of the strong coupling 
the least precisely known of all fundamental constants in nature. Improving our knowledge of a s 
is a prerequisite to reduce the uncertainties in perturbative-QCD calculations of all (partonic) 
cross sections at hadron colliders, and for precision fits of the Standard Model (a s dominates e.g. 
the Higgs boson H —> bb partial width uncertainty). It has also far-reaching implications includ¬ 
ing the stability of the electroweak vacuum 2 or the scale at which the interaction couplings unify. 

Having at hand new independent approaches to determine a s from the data, with experi¬ 
mental and theoretical uncertainties different from those of the other methods currently used, 
is crucial to reduce the overall uncertainty in the combined a s world-average value 1 . In Refs. 3,4 
we have presented a novel technique to extract a s from the energy evolution of the moments of 
the parton-to-hadron fragmentation functions (FF) computed at approximate next-to-leading- 
order (NLO*) accuracy including next-to-next-to-leading-log (NNLL) resummation corrections. 
This approach is extended here to include full NLO, plus a set of NNLO, corrections. Our 
new NNLO*±NNLL theoretical results for the energy dependence of the hadron multiplicity 
and mean value of the FF are compared to jet fragmentation measurements in e + e~ and deep- 
inelastic e^p collisions, in order to extract a high-precision value of a s . 

2 Evolution of the parton-to-hadron fragmentation functions at NNLO*±NNLL 

The distribution of hadrons in a jet is encoded in a fragmentation function, D^^z, Q ), which 
describes the probability that the parton i fragments into a hadron h carrying a fraction 
z = Phadron/Pparton of the parent parton’s momentum. Usually one writes the FF as a function 
of the log of the inverse of z, £ = ln(l/z), to emphasize the region of relatively low momenta 
that dominates the jet hadronic fragments. Starting with a parton at a given energy Q, its 




evolution to another energy scale Q' is driven by a branching process of parton radiation and 
splitting, resulting in a jet “shower”, which can be computed perturbatively using the DGLAP 5 
equations at large z > 0.1, and the Modified Leading Logarithmic Approximation (MLLA) 6 , 
resumming soft and collinear singularities at small z. Due to colour coherence and gluon radi¬ 
ation interference, not the softest partons but those with intermediate energies multiply most 
effectively in QCD cascades, leading to a final FF with a typical “hump-backed plateau” (HBP) 
shape as a function of £ (Fig. 1), which can be expressed in terms of a distorted Gaussian: 

£>(£, Y, A) = Af/(aV2^) • e [§ fe -§ s M( 2 + fc B 2 + 4 s.5 3 +ifcS 4 ] ^ with 6 = ^ (i) 


where A f is the average hadron multiplicity inside a jet, and £, a, s, and k are respectively the 
mean peak position, dispersion, skewness, and kurtosis of the distribution. The set of integro- 
differential equations for the FF evolution combining hard (DGLAP, MLLA, next-to-MLLA) and 
soft (DLA) radiation can be solved by expressing the Mellin-transformed hadron distribution in 


for t = In Q, leading to a 

0(aj )+0(a s )+0(aJ )+0(a s )+0(aJ )+•••• 


terms of the anomalous dimension 7 : D ~ C(a s (t)) exp f t 'y(a s (t , ))dt 
perturbative expansion in half powers of a s : 7 

Corrections up to order aj 2 were computed for the first time in Ref. 3,4 . The full set of NLO 
O(a^) terms for the anomalous dimension, including the two-loop splitting functions, plus a 
fraction of the 0(a s ) terms, coming from the NNLO expression for the ct s running, have 
now been computed 7 . Upon inverse-Mellin transformation, one obtains the energy evolution 
of the FF, and its associated moments, as a function of Y = ln(F//A QCD ), for an initial parton 
energy E, down to a shower cut-off scale A = ln((5 0 /A qCD ) for Nf = 3,4,5 quark flavors. 
The resulting formulae for the energy evolution of the moments depend on A QCD as single free 
parameter. Particularly simple expressions are obtained in the limiting-spectrum case (A = 0, 
i.e. evolving the FF down to Qo = A QCD ) motivated by the “local parton hadron duality” 
hypothesis for infrared-safe observables which states that the distribution of partons in jets are 
simply renormalized in the hadronization process without changing their HBP shape. Thus, by 
fitting the experimental hadron jet data at various energies to Eq. (1), one can determine cc s 
from the corresponding energy-dependence of its FF moments. 


3 Data-theory comparison and a s extraction 

The first step of our procedure is to fit to Eq. (1) all existing jet FF data measured in e + e _ and 
e , zz-p collisions at y/s ~ 1-200 GeV (Fig. 1) in order to obtain the corresponding FF moments 
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Figure 1 - Charged-hadron distributions in jets as a function of £ = ln(l/z) measured in e + e _ at y^s « 2-200 GeV 
(left) and e ± ,z'-p (Breit frame, scaled up by x2 to account for the full hemisphere) at -^s « 4-180 GeV (right), 
individually fitted to Eq. (1), with the hadron mass corrections (m e g = 130,110 MeV) quoted. 





















at each jet energy. Finite hadron-mass effects in the DG fit are accounted for through a rescaling 
of the theoretical (massless) parton momenta with an effective mass m e g as discussed in Refs. 3,4 . 
The overall normalization of the HBP spectrum (/C cll ), which determines the average charged- 
hadron multiplicity of the jet, is an extra free parameter in the DG fit which, nonetheless, plays 
no role in the final A QCD value given that its extraction just depends on the evolution of the 
multiplicity, and not on its absolute value at any given energy. 




Figure 2 - Energy evolution of the jet charged-hadron multiplicity (left) and FF peak position (right). Top: 
Comparison of theoretical predictions at four levels of accuracy. Bottom: Fit of the experimental e + e“ and DIS 
jet data to the NNLO*+NNLL predictions. The obtained IC c h normalization constant, the individual NNLO* 
tts(m^) values, and the goodness-of-fit per degree-of-freedom x 2 /ndf, are quoted. 


Once the FF moments have been obtained, one can perform a combined fit of them as 
a function of the original parton energy (which in the case of e + e _ corresponds to half the 
centre-of-mass energy y/s/2 and, for DIS, to the invariant four-momentum transfer Q DIS ). The 
top panels of Fig. 2 show the energy evolution of the zeroth (multiplicity) and first (peak 
position, closely connected to the mean of the distribution) moments of the FF, computed at an 
increasingly higher level of accuracy (from LO to NNLO*). The hadron multiplicity and FF peak 
increase exponentially and logarithmically with energy, and the theoretical convergence of their 
evolutions are very robust as indicated by the small changes introduced by incorporating higher- 
order terms. The corresponding data-theory comparisons are seen in the bottom panels of Fig. 2. 
The NNLO*+NNLL limiting-spectrum (A = 0) predictions for Nf = 5 active quark flavour^, 
leaving A QCD as a free parameter, reproduce very well the data. The most “robust” FF moment 
for the determination of A QCD is the peak position £ max which is quite insensitive to most of the 

“The moments of the lowest--^ data have a few-percent correction applied to account for the (slightly) different 
(Nf = 3,4) evolutions below the charm and bottom production thresholds. 










uncertainties associated with the extraction method (DG and energy evolution fits, finite-mass 
corrections) 4 as well as to higher-order corrections. The hadron multiplicities measured in DIS 
jets appear somewhat smaller (especially at high energy) than those measured in e + e _ collisions, 
due to limitations in the FF measurement only in half (current Breit) e^p hemisphere and/or 
in the determination of the relevant Q scale 4 . The value of a a (m^) obtained from the combined 
multiplicity+peak fit yields a s (m^) = 0.1205T0.0010, where the error includes all sources of 
uncertainty discussed in Ref. 4 . An extra theoretical scale uncertainty of Iq'qqqq should be added 
(this is a conservative estimate obtained only at NLO* accuracy 4 ). In Fig. 3 we compare our 
extracted a s (m^) value to all other NNLO results from the latest PDG compilation 1 , plus that 
obtained from the pion decay factor 8 , and the top-quark pair production cross sections 9 . The 
precision of our result (^S) is clearly competitive with the other measurements, with a totally 
different set of experimental and theoretical uncertainties. A simple weighted average of all these 
NNLO values yields: a s (m^) = 0.1186 ± 0.0004, in perfect agreement with the central value of 
the current world-average, a s (m^) = 0.1185T0.0006, but with a 30% smaller uncertainty. 
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Figure 3 - Summary of q s determinations using different methods at NNLO^ accuracy. The dashed line and 
shaded (yellow) band indicate the current world-average and uncertainty (listed also on the top) 1 . 
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